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We study the dynamics of a harmonically trapped quasi-one-dimensional Bose-Einstein condensate 
subjected to a moving disorder potential of finite extent. We show that, due to the inhomogeneity of 
the sample, only a percentage of the atoms is localized at supersonic velocities of a random potential. 
We find that this percentage can be sensitively increased by introducing suitable correlations in the 
disorder potential such as those provided by random dimers. 
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I. INTRODUCTION 

Suppression of wave transport in non-dissipative lin- 
ear systems can be induced by the presence of disorder: 
the scattered waves from the modulation of the random 
disorder potential destructively interfere in the forward 
direction, with a resulting vanishing wave transmission. 
This phenomenon is called Anderson localization (AL) 
[HQ. I n three dimensions (3D), AL takes place for states 
with energy less than a threshold (mobility edge). In 
two (2D) and one dimensions (ID) and in the absence 
of interactions, all single-particle quantum states are ex- 
pected to be localized Q. However, in the presence of 
correlations, the situation differs and a subset of delocal- 
ized states can appear in the spectrum |4| or an effective 
or even a true [l(| mobility edge can be observed 
even at low dimensions. 

Anderson localization of noninteracting atomic matter 
waves was observed in momentum and real space. In mo- 
mentum space using kicked-rotor setups in ID 11| and 
3D [13], while in real space in ID [rJ [l4j], and very re- 
cently in 3D [H [3. On the other hand, in 2D only 
anomalous diffusion has been observed [ItJ ■ In the last 
few years it has been experimentally demonstrated that 
disorder strongly affects the transport properties and dy- 
namics of a BEC, as for instance illustrated in (l8| . 

One of the outstanding challenges of physics is to un- 
derstand the interplay between disorder and interactions. 
In the case of an interacting condensate, wave scattering 
from the random potential does not occur if the wave 
group velocity is lower than a critical velocity v c that co- 
incides with the (local) sound velocity c in the limit of 
small random potential amplitude and decreases down to 
vanishing values in the strong disorder limit [l9j . Thus, 
to observe AL in an interacting BEC it is necessary that 
the speed of the relative motion between the superfiuid 
and the disorder potential is larger than v c . This s etup 
was proposed and theoretically studied in Refs. |20l422| . 
These authors studied the flow of an homogeneous quasi- 
one-dimensional Bose-Einstein condensate through a dis- 



order potential of finite extent. That disorder potential 
moves with a velocity v with respect to the condensate. 
In the subsonic regime, the flow is superfiuid and the 
density profile is stationary. In the opposite supersonic 
regime, a region of stationary flow also exists, but in this 
case energy dissipation occurs. In this domain, depend- 
ing on the extent of the disorder potential, the system is 
either in a Ohmic or in an AL regime, respectively char- 
acterized by a transmission decreasing linearly or expo- 
nentially with the size of the system L. 



At variance of Refs. [20 22] , in this paper we study the 
effects of the inhomogeneity of a cigar-shaped trapped 
BEC in the presence of a moving disorder potential. We 
investigate the possibility of observing BEC localization 
by looking at the position of the center of mass of the 
condensate. If the center of mass moves along with the 
moving potential then the system shall be in the AL 
regime or in another kind of localized phase. Because 
of the inhomogeneity, we observe the localization of only 
a percentage of the atoms in the BEC. This percentage 
of localized atoms can be increased or suppressed by in- 
troducing ad-hoc short-range correlations on the random 
potential. 

The paper is organized as following. In Sec. [TT] we 
introduce the time-dependent nonpolynomial nonlinear 
Schrodinger equation (NPSE) that describes the conden- 
sate dynamics in the elongated geometry and in the pres- 
ence of a moving disorder potential that we characterize 
by its auto-correlation function. As discussed and shown 
in Sec. IIII1 the disorder potential drags the atoms with an 
efficiency that depends on both the scattering properties 
of each impurity and on the impurity density. The case 
of two types of impurities, single and dimerized, at dif- 
ferent densities, have been studied, highlighting the role 
of correlations in the localization dynamics. Because of 
the inhomogeneity of the BEC, we observe that the drag 
force is more efficient in the BEC tails where the local 
sound velocity is lower and superfluidity breaks down at 
small drift velocities v. Our concluding remarks are given 
in Sec. HV1 



II. THE MODEL 
A. Equation of motion for the BEC wavefunction 

Our starting point is the equation of motion for a 3D BEC 
trapped into a cigar-shaped potential. Such an equation 
is known as the Gross-Pitaevskii equation (GPE) 
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The wave function ip(r,t) describes the condensate, 
which is constituted of N atoms of mass m. g — 
Anh 2 a s /m stands for the interaction coupling constant 
with a s being the s-wave scattering length between 
atoms; in our system, interatomic interactions are repul- 
sive and, as a result, a s > 0. The trapping potential U(r) 
is given by the sum of a static cigar-shaped harmonic trap 
and a time-dependent random potential: 
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U(r, t) = -muji(x 2 + y 2 ) + -muj 2 z z 2 + V(z, t) (2) 

with wj_ and ui z the trapping frequencies in the per- 
pendicular and longitudinal directions, respectively and 
lj z <C The last time-dependent term in ^ corre- 
sponds to a random potential that is fixed in the moving 
frame z — z — vt, v = ve z being the drift velocity. 

Under this trap geometry and a further assumption 
discussed below, the 3D GPE can be reduced to an ef- 
fective ID time-dependent NPSE [23|. The advantage of 
the ID NPSE is that it is easier to deal with when mak- 
ing computations. In order to obtain such dynamical 
equation, we begin with a variational Ansatz 

-(x 2 +y 2 )/2a 2 (z,t) 

^(r, t) = f(z, t)cf>(r, t) = f(z, t) — — ^ (3) 

where the transverse part <p(r, t) is modeled by a Gaus- 
sian function with variance a(z,t). The validity of this 
description is based on the assumption that a(z, t) slowly 
varies as a function of z and t such that the kinetic energy 
term d 2 /dz 2 associated to ^(r, t) can be neglected. Both 
the longitudinal wave function f(z,t) and the variance 
a(z,t) are determined by the energy variational princi- 
ple. For f(z,t) one gets the NPSE 



<4/ 

dV 



& d 2 1 22 „. . 
'2^+2™^ +V{2 ' t) 



-tvjj i 



l + 3a s A|/| 2 
^1 + 2aJV|jf 



(4) 



The variance is given by 



a 2 (z,t) = ai^l + 2a s N\f{z,t)\ 2 , 



(5) 



where a± = yjh/ (moj±) is the oscillator length in the 
transverse direction. The 3D density profile and velocity 




FIG. 1: Schematic representation of the disorder potential. 

(a) The Gaussian peaks are randomly distributed (AM). 

(b) The peaks are distributed in a random dimer sequence 
(RDM). 



field are then 
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l s « f*(z)f(z)-f'(z)r(z) ^ 
v(r)=w(r)z = — z, (7) 
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with p(z) = \f\ 2 the integrated ID density. 

The NPSE is numerically solved using a split-step 
method and spatial Fast Fourier transforms (FFT). First 
we compute the equilibrium density profile in the pres- 
ence of a static disorder potential. Then, we switch on 
the drift velocity v and compute the time evolution of 
the condensate wavefunction f(z,t). In this work, we fo- 
cus on a system of 10 5 87-Rubidium atoms subject to a 
transverse confinement of oj± — 2tt x 500 Hz and a lon- 
gitudinal confinement of lo z = 2n x 7 Hz. The s-wave 
scattering length has been fixed at a s = 80 Bohr radii. 



B. The random potential 

The random potential V(z, t) is modeled by the sum 
of -/Vdis Gaussian functions of height Vdi s and width w, 
randomly distributed at positions zi — jid, where ji is a 
random integer number and d fixes the minimal distance 
between the peaks. Such a disorder potential could be 
realized by deeply trapping some impurities (heavy atoms 
of another species) in an optical lattice strongly detuned 
from the condensate atomic frequencies 24-26]. When 
the disorder pattern is pulled with a constant speed v 
through the system, if v is lower than the sound speed 
c = a/m7 (2m) (23], we expect the disorder not to affect 
the system because of the superfluid nature of the gas 
itself [ifj. In contrast, the effect of the disorder potential 
should appear at v > c where the kinetic energy starts 
to compete with the interaction energy and the limit of 
a noninteracting gas is reached for «»c. 

We will consider two sorts of potential patterns: (i) 
an Anderson-like distribution, that we will call Anderson 
Model (AM), where the jVs are randomly distributed; (ii) 
a Random Dimer Model (RDM) distribution where peaks 
are dimerized and the dimers are randomly distributed 
(see Fig [I}. The dimers are stamped by the peak-to-peak 
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FIG. 2: (Color online) Rescaled auto-correlation function 
C(\zi - z 2 \)/C(0) as a function of (zi - z 2 )/X for the AM 
(red continuous line) and RDM (green dashed line) for the 
case d = A. The top panel corresponds to ridis = 0.12A - , 
and the bottom panel to ridis = 0.5A -1 . 
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FIG. 3: (Color online) Same as in Fig. [2] but for the case 
d = A/2 and ridis = 0.5A -1 . 



dimer structure with £ — X. This is well visible at low 
peak density (top panel), while by increasing the bump 
density the probability to find dimerized structures in 
the AM potential increases as well. This is seen in the 
heights of the peaks at z\ — z 2 — jX for the two models 
becoming closer even for j = 1 (bottom panel). 

In Fig. [3] we show C(\zi — Z2I) for the case d = A/2 
and ridis = 0.5A -1 (keeping fixed £ = A). In this case 
the RDM peak at z\ — Z2 = X is lower of more than 
a factor 2 with respect to the case with d = X at the 
same disorder density (bottom panel of Fig. [5]) : we thus 
expect the dimer structure to play a minor role for the 
case d = A/2. 



distance £. Hereafter we consider a disorder potential 
characterized by an amplitude Vdi s = 0.02E r , where the 
recoil energy E r — h 2 / (2mA 2 ) for each atom of mass m is 
defined with respect to the wavelength A = 780 nm, char- 
acterizing the D2 hyperfine Rubidium transition. The 
dimer peak-to-peak distance £ was set equal to A and the 
width of a single bump w is fixed at 140nm, roughly A/5, 
ensuring no sizeable overlap between Gaussian functions. 
The disorder potentials can be characterized in term of 
their autocorrelation functions 

C(\ Zl -z 2 \) = ((V(z 1 ,t)-(V(z 1 ,t)))(V(z 2 ,t)-(V(z 2 ,t)))), 

(8) 

which strongly influence the nature of the energy states 
for low amplitude disorder near equilibrium. In the case 
of the AM and RDM potentials in Fig. fj we plot the 
autocorrelation function for the case d = X and for two 
different peaks densities: ndi s = 0.12A -1 (top panel) and 
ridis = 0.5A -1 (bottom panel). It is worthwhile remark- 
ing than the peak density corresponds to the average 
number of Gaussian peaks both in the AM and RDM 
potential and thus in average the number of dimers in 
RDM is half the number of peaks in the AM. The mod- 
ulation with spatial period A for both the AM and the 
RDM evinces that bumps and dimers are both randomly 
distributed over discretized positions of step d — X. The 
main difference between the AM and the RDM is that 
the RDM has a larger peak at z\ — z 2 = A because of the 



III. THE LOCALIZATION FRACTION AND 
THE DRAG FORCE 

If the condensate or a part of it is localized, we expect 
it to follow the pulled disorder potential. The localized 
portion of the condensate glues to the disorder "band- 
wagon" and therefore travels the same distance as the 
disorder potential. This dynamics depends on the forces 
experienced by the atoms. The force acting on the BEC 
center-of-mass has two terms, F = + i<di s , one due 
to the harmonic confinement, = —muj 2 z c . m . and the 
other due to the disorder potential 

r + 00 fly 

F ^ = -J dz\f(z,t)\ 2 — . (9) 

For small center-of-mass displacements Az c . m . ~ 0, the 
leading term is the drag force F^- ls due to the disorder 
potential. In this regime the localization fraction Ni oc /N 
can be deduced by the ratio between the Az c m . and the 
distance Azdi s traveled by the disorder potential in the 
same time interval, namely in this regime we can identify 

Nioc _ Az c . m , 
N ~ A^is ' 1 ' 

Indeed if the center of mass travels the same distance as 
the disorder potential, this would mean that the 100% 
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FIG. 4: (Color online) Localized BEC fraction (in logarithmic 
scale) as a function of the disorder potential drift velocity v, 
in units of the sound speed c evaluated at the center of the 
trap. The red crosses correspond to a single bump and the 
green stars to a dimer. 



of the atoms are localized. The localized condensate will 
stop following the moving potential and it will change 
direction at a time tj at the position Zf verifying 



Fdi s (t)v c . m .{t)dt 
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Thus, the turning point Zf will provide a direct measure 
of the average value of the drag force during the conden- 
sate forward motion via the relation 

z f = 2F dis /(rmj 2 z ), (12) 

where we have defined 



Fa 



F dis (t)v c . m .(t)dt/zf 



F dis (t)dz c , m ./z f . 

(13) 

The dependence of the drag force with the drift velocity 
gives us a direct measure of the loss of superfluidity in 
the system. According to the Landau criterion |28| a 
single impurity is expected to flow without friction below 
a certain velocity, corresponding to the sound velocity for 
a weakly interacting BEC. 



A. The single impurities 

The localization efficiency of a disorder pattern de- 
pends on the impurity density and on the reflectivity of 
each impurity. In this work we are comparing single im- 
purities randomly distributed with dimerized structures. 
With the aim to understand the difference in behavior 
of the localization efficiency of a dimer with respect to a 
single bump, we first look at the condensate dynamics in 
the presence of a moving single bump (red crosses in Fig. 
2]) and of a moving dimer (green stars in Fig. 2]) . 

In Fig. @] we plot the fraction of atoms that follows 
the moving defect over a distance of 133/im. We observe 




FIG. 5: (Color online) Analogous optical system of the dimer 
structure present in our model. 



that the localization fraction of the dimer is a strongly 
non-monotonic function of v. The single dimer localizes 
the atoms less efficiently than a single bump for v ~ 0.9c, 
while it is more efficient by a factor of 3 over a velocity 
range of 1-1.2 v/c. The suppression and the enhance- 
ment of the localization are both a signature of some 
interference effect due to the internal structure of the 
single dimer. This behavior can be qualitatively repro- 
duced by considering the analogous optical system of our 
model, schematically shown in Fig. [5] It consists of 
two dielectric slabs of refraction index n' , width w, at 
distance a — I — w, merged in a medium of refractive 
index n, with n = 1, and n' ~ 1 + V d i S /mv 2 in the limit 
ran 2 3> Vdi S (and v > c). This model corresponds to 
associate to an incident wave of energy E = mv 2 /2 and 
wavevector k = mv/h a transmitted wave of wavevec- 
tor k' = Ky/l — V d i S /E ~ k(1 — Vdis/2E) in the regions 
where the disorder potential is present. The reflection 
coefficient for an incident wave of wavevector k through 
the two-slab system (from region 1 to region 5, as shown 
in Fig. [5]) can be written as 



with 
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where a = ^?kw, f3 = na, r 34 = r 12 = -r 23 = -r 45 = 
(n — n')/(n + n 1 ). Equation (IT4l) takes the form 

ri5 = r i2 x 

[-1 - 2e 2M + e 2 ^ + rf 2 e 4fo + r 2 2 e 2l <- a +^ - 2e 2l ^+ 2Q )] 
1 - 2rl 2 e 2ia + r 2 2 e 4M + r 2 2 e 2 ^ - 2r 2 2 e 2l ( Q +« + 2rf 2 e 2i ( 2a +^ 

(16) 

The behavior of the reflectivity 1Z = |ris| 2 of the dimer 
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FIG. 6: (Color online) Reflectivity (in logarithmic scale) of a 
single defect (red continuous line) and of a dimerized structure 
(green dashed line) as a function of the disorder potential 
velocity. 

structure must be compared with that of a single bump 



This is shown in Fig. |U For our choice of the parame- 
ters, the reflectivity of the dimer oscillates with respect 
to that of a single bump that decreases monotonically, 
in qualitative agreement with what observed for the lo- 
calization fraction shown in Fig. HI The shift of the 
minimum position for the dimer reflectivity with respect 
to the localization may be attributed to on the one hand, 
the Gaussian shape of bumps as compared to the rect- 
angular shape of the dielectric slabs, and on the other 
hand to the inhomogeneity and non-linearity of the sys- 
tem. All these factors are not taken into account in the 
current optical model. Finally, let us remark that both 
the single bump and the single dimer yield full dereal- 
ization (1Z = 0) for a = 7r, namely when each bump plays 
the role of a cavity [1, [26| , corresponding to v c in our 
system. 

B. Random distribution of impurities 

As already illustrated in Sec. Ill Bt we consider three 
sets of parameters for the AM and RDM potentials: (i) a 
bump density rtdis — 0.12A -1 ~ 0.16 pcaks/^tm, and d = 
A; (ii) a bump density rtdis — 0.5A -1 ~ 0.65 peaks/^m, 
and d = A; and (iii) with the same bump density as in 
(ii) but with d = A/2. In the three cases, the size of the 
disorder potential is 485 /mi and the BEC size is 204 //m. 

We run the simulations for up to approximately 3-4 
cycles of ui z . We generally note that in our range of 
drift velocities, longer time durations are not necessary 
since the moving disorder potential would try to bring 
the condensate too high up along the harmonic potential 
and thus the condensate would inevitably fall back at Zf. 

We run simulations for 37 different drift velocities, in 
the range v = [0.46 - 1.30]c = [1.40 - 3.98]mm/s (where 



c = 3.06mm/s), with between 3 and 24 random poten- 
tial realizations for each drift velocity. To obtain the 
percentage of localized atoms, we compute the distance 
traveled by the condensate center-of-mass correspond- 
ing to a fixed value of the potential drift. We used 
that unique velocity- independent distance, in our case 
Azdis = 96/im, in order to compare equivalent final po- 
tential configurations. For this value of Azdi s , we mea- 
sure a small center-of-mass displacement Az c m ., gener- 
ally lower than an oscillator harmonic length in the axial 
direction, Az c m . < a z — ^Jh/(muj z ). Then we identify 
the ratio Ni oc /N with the ratio between the center-of- 
mass shift and A^dis according to Eq. (fTUj) . Figure [Jj 
shows the localized ratio Ni oc /N for the potentials (i), 
(ii) and (iii). We first note that irrespective of the kind 
of disorder the localization occurs for v > 0.5c. This 
unexpected localization regime may be due to the local 
sound velocity being inhomogeneous due to the harmonic 
confinement [29j or to the non- vanishing amplitude of the 
disorder potential [19j . Ni oc /N is very small for all data 
sets, except for the parameters (ii) at v ~ 1.1c, for both 
AM and RDM (see the middle plot in Fig. [7]). We claim 
that this is an effect related to the presence of dimers 
both in the RDM by construction, and in the AM due 
to the choice d — A and density ridis — 0.5A -1 , as shown 
in the bottom panel of Fig. [2j Indeed, in the low den- 
sity case (top panel of Fig. [JJ where the AM has just 
developed few dimers with respect to the RDM, the lo- 
calization enhancement in the RDM is well visible, even 
if the localization fraction is very small. Moreover, its 
oscillatory behavior as function of v reminds that of the 
single dimer shown in Fig. [4j with minimal localization 
near v/c ~ 0.9. Finally, in the bottom panel of Fig. [JJwe 
show that the choice d = A/2 suppresses the Ni oc /N peak 
at v ~ 1.1c, since it destroys the correlations introduced 
by the single dimer itself (cf. Fig. [3]). 

In Fig. [5] we plot the turning point Zf for the AM 
with ndis = 0.5A -1 and d — A/2 and the RDM with 
the same peak density but with d = A. We verified nu- 
merically that for the disorder amplitude and velocities 
considered the turning point Zf is proportional to the av- 
erage drag force Fdis as expected from Eq. (fT2"]l in the 
limit of small displacements. Not surprisingly, the over- 
all behavior of the average drag force Fdi s is qualitatively 
similar to the localization ratio, showing that to a greater 
drag force, i.e., less superfluid fraction, it corresponds a 
higher localization efficiency. In agreement with the sin- 
gle defect analysis (Sec. MI A[) . in both Figs. [7J and [51 we 
clearly observe a partial suppression of the localization 
at v ~ 0.8 — 0.9c, and an enhancement of the localiza- 
tion at a supersonic region of the drift velocities, when 
correlations arc dominant. 



C. Localization inhomogeneity 

In order to prove that we are observing the localiza- 
tion of a BEC fraction rather than a slower drag of the 
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FIG. 7: (Color online) Localized BEC fraction as a function 
of the disorder potential drift velocity v, in units of the sound 
speed c evaluated at the center of the trap. The red crosses 
correspond to the AM and the green stars to the RDM. The 
top panel corresponds to rials ~ 0.12A -1 , and d = A; the 
middle panel corresponds to ridis — 0.5A" 1 , and d — A; the 
bottom panel corresponds to n^is — 0.5A" 1 , and d = A/2. 



whole BEC, we analyze the dynamics of the BEC tails, 
where the density is lower and AL should occur more ef- 
ficiently. In particular, we focus on the forward moving 
tail distribution that experiences the disorder potential 
for the whole simulation time. Analogously to Sec. IIIIB1 
we compute the percentage of atoms localized in the for- 
ward moving tail of the condensate Ni oc j/N t by dividing 
the tail center-of-mass shift with the same Azdi S value. 
We looked at 3 different sizes of that tail: those compris- 
ing 1%, 5% and 10% of the total condensate mass. In 
the case of the forward tails that amounts to 1% of the 
total condensate mass (see Fig. [S]) , we observe an almost 
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FIG. 8: (Color online) Turning point Zf as a function of the 
disorder potential velocity v, in units of the sound speed c 
evaluated at the center of the trap. The red crosses corre- 
spond to the uncorrelated disorder with ridis = 0.5A -1 and 
d — A/2, and the green stars to the correlated one with d = A 
and the same peak density. 
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FIG. 9: (Color online) Localized BEC fraction of the leading 
moving tail -1% of the total condensate mass- as a function of 
the disorder potential velocity v, in units of the sound speed 
c evaluated at the center of the trap. The red crosses corre- 
spond to the AM (hi) and the green stars to the RDM (ii). 



complete localization in that same supersonic region of 
the drift velocities. While for the 5% and 10% mass tails 
(not shown), we measure 50% of localization efficiency. 
The negative values of localization in the subsonic region 
of drift velocities are due to the fluctuations of the con- 
densate density caused by the presence of the disorder 
potential. At those velocities, these fluctuations provoke 
the forward leading small portion of mass to buckle back 
towards the center of the condensate even though there 
is no much overall motion of the whole condensate. 



IV. CONCLUSIONS 

We have presented an analysis of the drag properties 
of two kinds of disorder potentials of finite extent mov- 
ing through an inhomogeneous quasi- ID Bose-Einstein 
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condensate. Because of the presence of the external 
harmonic confinement our system, unlike (2lj . is never 
in a stationary state. We treated both cases of non- 
correlated and correlated disorder with short-range cor- 
relations. Our numerical computation of the fraction of 
localized atoms and the drag force shows that for the 
case of correlated disorder, in the form of random dimers, 
there is a suppression or an enhancement of localization 
depending on the drift velocity. This was buttressed by 
our analytical optical model in which we determined the 
reflectivity of a single dimer and of a lonely defect. The 
effects of correlations arc masked as we increase the disor- 



der density and the dimers cannot be distinguished from 
a high density collection of single defects. 
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